{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Heat Capacity of Cementite ($Fe_3C$)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Bengt Hallstedt, Dejan Djurovic, Jörg von Appen, Richard Dronskowski, Alexey Dick, Fritz Körmann, Tilmann Hickel, Jörg Neugebauer, Thermodynamic properties of cementite, Calphad, Volume 34, Issue 1, March 2010, Pages 129-133, ISSN 0364-5916, http://dx.doi.org/10.1016/j.calphad.2010.01.004. (http://www.sciencedirect.com/science/article/pii/S0364591610000052)\n",
    "\n",
    "The TDB file used here differs slightly from the published TDB to ensure compatibility with pycalphad's TDB parser. All phases except cementite are omitted. The numerical results should be the same."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "TDB = \"\"\"\n",
    " ELEMENT C    GRAPHITE                   12.011     1054.0      5.7423 ! \n",
    " ELEMENT FE   BCC_A2                     55.847     4489.0     27.2797 ! \n",
    " TYPE_DEFINITION % SEQ * !\n",
    " TYPE_DEFINITION A GES AMEND_PHASE_DESCRIPTION @ MAGNETIC -3 0.28 !\n",
    " PHASE CEMENTITE_D011 %A 2 3 1 !\n",
    " CONSTITUENT CEMENTITE_D011 : FE : C : !\n",
    " PARAMETER  G(CEMENTITE_D011,FE:C;0) 0.01 +GFECEM; 6000 N !\n",
    " PARAMETER  TC(CEMENTITE_D011,FE:C;0) 0.01 485.00; 6000 N !\n",
    " PARAMETER  BMAGN(CEMENTITE_D011,FE:C;0) 0.01 1.008; 6000 N !\n",
    " FUNCTION GFECEM      0.01  +11369.937746-5.641259263*T-8.333E-6*T**4;\n",
    "        43.00  Y  +11622.647246-59.537709263*T+15.74232*T*LN(T)\n",
    "       -0.27565*T**2;\n",
    "       163.00  Y  -10195.860754+690.949887637*T-118.47637*T*LN(T)\n",
    "                  -0.0007*T**2+590527*T**(-1);\n",
    "      6000.00  N !\n",
    "\"\"\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Do some initial setup, including reading the database."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Only needed in a Jupyter Notebook\n",
    "%matplotlib inline\n",
    "# Optional plot styling\n",
    "import matplotlib\n",
    "matplotlib.style.use('fivethirtyeight')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from pycalphad import Database, calculate\n",
    "\n",
    "db = Database(TDB)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Compute the molar heat capacity at all temperatures from 1K to 2000K with a step size of 0.5K.\n",
    "\n",
    "We do this with the `calculate` routine instead of `equilibrium` because the cementite phase has zero internal degrees of freedom. Since there's nothing to minimize, we can do the computation faster with `calculate`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "result = calculate(db, ['FE', 'C'], 'CEMENTITE_D011', T=(1, 2000, 0.5), output='heat_capacity')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAGTCAYAAAClNV9VAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlYVPX+B/D3mRn2RUAR9bqQgmuIpgai5tJNJVQytbTu\ntVzS1OtVuy5RZmaaZZqaEVZm2c3S9GKZlqXpbXHtVwZ6VcRUREWQZdiXYeb8/kAGzsCMB5hhhpn3\n63l68Cwz8zm48O67Cmq1WgQRERERNXkKaxdARERERObBYEdERERkJxjsiIiIiOwEgx0RERGRnWCw\nIyIiIrITDHZEREREdoLBjoiIiMhOMNhZQHJysrVLsCo+P5/fkfH5+fyOjM9v/ednsCMiIiKyEwx2\nRERERHaCwY6IiIjITjDYEREREdkJBjsiIiIiO8FgR0RERGQnGOyIiIiI7IRK7o1//PEHTpw4gaSk\nJGRlZUEQBDRv3hydO3dGWFgYevfubck6iYiIiOguTAa727dv44MPPsDnn3+OGzduQBRFODs7w8fH\nB6IoIjc3F2VlZRAEAW3atMGkSZPwzDPPoGXLlo1VPxERERHdYTTYvfLKK3j//ffh5eWFMWPGYMiQ\nIejduzdat24tuS8tLQ1//PEHDh8+jH//+9949913MWPGDLz88ssWL56IiIiIqhgNdj/++CM2b96M\nUaNGQRAEo2/QunVrtG7dGpGRkVizZg2+/vprbNiwwSLFEhEREZFxRoPd4cOH6/xmgiBgzJgxGDNm\nTIOKImoK0oq0ePtMPpwUAp7r6QUfF85FIiIi6zL5k+jSpUt1erMVK1Y0qBiipqJcJ2Lcd5mIO1eI\nt88WYN6xHGuXREREZDrYjR49Gn/++aesN/rXv/7FLlhyGDv+LMI5dbn++PCNUitWQ0REVMFksNPp\ndHcNdzqdDjNmzMDWrVvx2GOPmb1AIlsjiiLWJ+ZLzuVrRCtVQ0REVMVksNu3b5/JcFdWVoYnn3wS\nu3btwjPPPIPNmzdbrFAiW3HkZin+zNNauwwiIqIaTAa74OBgfP3117WGu8LCQowfPx4HDhzAwoUL\nsWbNGosXS2QL1hm01gGAm9L4zHEiIqLGctdpfMHBwdi3bx9EUdSHO7VajejoaPz8889YuXIlXnzx\nxcaolcjqTqaX4uitshrnRbArloiIrE/WlmJBQUHYt28fRo8ejVGjRsHHxwfJycl4++238fe//93S\nNRLZjNpa6wAw1hERkU2QvfBWp06dsH//figUCly5cgVbt25lqCOHcuxWKb6/XvvsVx2THRER2QCT\nLXatWrWqseuERqOBTqfDzJkzMXPmTMk1QRBw8+ZN81dJZGWiKGLFb3kmrjdiMUREREaYDHZjx441\nuZ0YkaM4kFqCExk1x9ZVYq4jIiJbYDLYxcXFNVYdRDarXFeztW54WxdJtyyDHRER2QJubkl0F1su\nFOJ8tV0mBADL+jST3MMxdkREZAvqHeyysrIQGhqKU6dOmbMeIpuSUazFa79LW+smBrmjh6+sCeVE\nRESNqt7BTqvV4tq1ayguLjZnPUQ25ZXf8pBXbbswbycBy/t41zr2VOQMCiIisjJ2xRIZcexWKbYn\nF0nOLentjQB3JYCKLtnqGOuIiMjaGOyIalFUrsM/fsmRnOvqo8KMbh76Y8NGO46zIyIia6t3sPPw\n8MCSJUsQGBhoxnKIbMNrv+fjcr5Wcu7NcB84KarSHFvsiIjI1tR7BLiHhweef/55c9ZCZBN+zSjD\nu+cKJOemdfXAoNYuknM1gh2THRERWVm9gl1+fj7y8vKg0+lqXGvXrl2DiyKylnyNDs/+nC3pVm3n\nqcTyvt417hUESJrpmOuIiMja6hTsPv74Y2zatAlXrlwxek92dnaDiyKylkXH1fgzT9oF+3aED7yc\nao5aUHCMHRER2RjZY+w++eQTLFiwAB06dMDSpUshiiJmzZqFBQsWoGXLlggJCcGmTZssWSuRRe24\nVIQdf0qX75naxQND/+Ja6/01x9gx2RERkXXJDnabN2/GkCFDEB8fj6effhoAMHz4cLz00ks4ceIE\n1Go18vKMb5JOZMsu5Wqw8Lhacq6bjwqr7m9m5BWAYBDtOMaOiIisTXawu3z5Mh5++OGKFykqXqbR\naAAAPj4+mDx5MrZs2WKBEoksK69Mhyd/yEZBeVUyc1UCHw7xg5uq5kLElQyXO2GuIyIia5Md7Dw8\nPPQr63t6ekKpVCItLU1/3c/PDzdv3jR/hUQWpBNFPPtzDpJyyyXnX7vfB919nUy+1vAvD4MdERFZ\nm+xg17lzZyQlJQEAVCoVQkJCsHPnTmg0GpSUlGDnzp3o0KGDxQolsoQ1f+Tjm2slknMTOrphShf3\nu76WCxQTEZGtkR3sHn74YRw4cAAlJRU/BBcuXIhjx44hMDAQQUFBOHnyJBYsWGCxQonMbe/VYrz+\nR77kXGhzJ7w9wLfWvWANcR07IiKyNbKXO5k7dy7mzp2rP46KisL+/fuxd+9eKJVKjBw5EgMHDrRI\nkUTmdiK9FDN+ki7N08JVgU+HmR5XJyHzNiIiosZS750nAKB///7o37+/uWohahTJuRpM+iELJdWW\nq1MJwMdD/dDOU/5fCcPm7k1n8/FSH+OzaImIyLJ0oohyHVAuitDoAK1ORLkIaHRAuU6EVqz4qrlz\nj1YHaO7cU/m68jv3losGX6tfF6veu/IerU7E/SoBwVb+HjQo2DXU0aNHsWnTJiQkJCAtLQ2xsbF4\n8sknAVTMuF25ciUOHjyIq1evwsvLC4MGDcLLL78s2d2itLQUS5cuxX/+8x+UlJTggQcewLp16/CX\nv/zFWo9FNiy9SItx32chp1Tab7o+wgcDW7kYeVXtuvg44WRGmf54XWIBlAoB4z3MUioRUYOJYkWY\nqT3I1Awr1c9VD0aVYUh7516NPsxU/7WIW7dVaFaYpz/WVHv/ijpE4wGqlhqMhqzKIGZQn7VHxLzR\nVfYIN4sxGuxGjx5d5zcTBAF79+6VfX9hYSG6d++OSZMm4dlnn5VcKyoqQkJCAhYuXIiQkBDk5eVh\n6dKlGD9+PI4ePQqVqqL0mJgYfPPNN/jwww/h6+uLF198EY8//jh+/PFHKJXKOj8D2S91qQ6PHcrC\ntQLpzhLP9/LC3zvXPY2t7e+DqG9uI09T9U/Jmj/y8XsLZ2zpoIOPi/X/ghNR3VSGEU21AKOpDC+S\nX9+5Vi1UaHQiUjOVSFQW3fU+rf68qG9NqvGayvOikfOV71VbK1S1MNW4nIGU/LvfZqca//tdk9Fg\np9PpZA0gr06s4+jx4cOHY/jw4QCA2bNnS641a9YMX375peTc+vXrER4ejqSkJPTo0QO5ubn497//\njdjYWAwdOhQA8N577yEkJAT//e9/8eCDD9apHrJfBRodHjuYhYQsjeT834LdsaSXV73eM8TPCZ/9\ntTnGf58p6dY9lKlCxJfpWNGvGcbd41bnv0dETZl4J2SU6USUaSuCR5k+kFQ/d+e8tiKQlFVe1wFl\n2qrQUnYnsFScq/5etYccya9Fw2BWy2uqtQqZp8XHBbiQY4bvJDVFNh3s9u/f35h1yJKfX/F/AT4+\nPgCAP/74AxqNBsOGDdPf07ZtW3Tp0gUnT55ksCMAQHG5iEmHsnDqdpnk/IN/ccH6CJ8GBa+BrVzw\n1YgWmPiDtHv3ZpEO03/MQez/CjC3hydGB7rByXBzWaJ6EMWKwFKiFVGmFVGiBUq1ouQ//Tld1bnK\ncKQPSjoRGm1VoKo8VyYJWnfCVbVzVWGsKmRJ38sN4lGuaUrW46QAVIIAlQJQCoCTovLXQsU1hQCV\nACgVd46FqnOSr9Xep+p8zXPKaq8JEtOt/fjWHWNXF2VlZVi6dClGjhypHz+XkZEBpVKJ5s2bS+71\n9/dHRkaG0fdKTk62aK2N9Rm2zFaev0wHLDrvgmM50m75nl5aLGufg6t/Nvz/rP0AbLlXQMwFF1ws\nlHa/ns7UYOqPOWhxPAsPNtfiry20uNdbB7kTb5sqW/n9N7fKsUqlOkAjAqU6oaKFSQeU6QSUiUCZ\nToEfT1yuOi/euaY/rvp1qU6oCEzinffUCTXeW1N53517ynSAaNNTsm25NsekhAilUBFyVHe+KgVA\npRChxJ1ziopAUHFNrAhCqLwPd+4TJcf69xDuXFdU+wzUfk0piFBJjitrEiX1VX01UrsgGhxXfFWg\n5hqjjc3S//4FB5ueniE72KWmpsq6r/rEBnMpLy/HjBkzkJubi88///yu94uiaLIV5m7flIZKTk62\n+GfYMlt5fo1OxNT/ZuNYjnQB4l7NnfDVyNZo5my+MXDBAH7uIWL16TxsOpsPrSj985dZpsDONAV2\npjnBy0lAeEtn9G/lgnt9ndDDzwlt3BV202Vrid9/URRRqoWkBaqyZaqi1aqi1ahEe+e+ynPaaud0\nIkrLq1qxJK81vE9yvuq+chvoZiHLUggVLT5Od1plKlt7VAoBTtVaf5wUVa0/1X9dWlQAXy+vatdq\nf42qlvdyqna+9s8xvEdaX60tTQIa9d8WW/n331ps4fllB7uePXvK+sORnZ1913vqory8HNOmTcO5\nc+ewb98++Pn56a+1bNkSWq0WWVlZaNGihf58ZmYmIiIizFoHNS2l2opQt99gV4nuPirED29u1lBX\nyUUpYHnfZhjkfBtbb/vW+OxK+RoRB2+U4uCNUv25Zs4COniq0M5TiXaeSrT1UKK5qxJ+Lgr4ugjw\nc1HAx0UBd5UAV6UARSP8Q60TK8ZDld7piivVVnS9lep/XTUeqvq11FtK+GgLUFpL6CrT4k6Aqvaf\nriqIVZ2v1r2oqzimpkMpAM4KAU7Kiq/Od4JH5a+dlHe+KgQ4K+8EJqWR65UBptp1laLqem3hRyXU\nEoYqfy3UEqaqHTf071ZycjaCg/3ufiORhcgOdu+8806NYKfVanHt2jXs2LED/v7+mD59ulmL02g0\nmDp1Ks6fP499+/YhICBAcr1Xr15wcnLCkSNHMGHCBADAjRs3kJSUhLCwMLPWQk1HcbmIyYezJMEJ\nAIK8VfhyZAv4uVp2tnR7NxHbH2yOJLUGm88V4Is/i1F4l6ae3DIRidkaJGZrTN5XyUUJuCoFuCkF\nuKoEqAQBCqGitaGyK0JR7ZxOrOhG1FabQaetPCdWnqtayqC0Qa1TLsCl3Pq+mGRQChW//853/hy4\nKAW4KO58VeLO14r/XO+EJJc7wUilqAxdFecNf1153VlZEXSc77QeOSurwo9zZeCqFtyclAKcBAFX\nL19C186O22JDZG2yg13l+nK1mT9/Pv76178iLy+vTh9eUFCAy5cvA6iYhXv9+nUkJibC19cXrVu3\nxlNPPYXTp0/j888/hyAISE+vGJTo7e0NNzc3NGvWDH//+9+xbNky+Pv765c76dGjB4YMGVKnWsg+\nFGp0mPRDNn5Kk4a6Dp5KfDWyBVq6Nd4SOF18nLA+whev3e+Dg9dLsDelGP+9WYrMEl2D37uyRSvX\n6qs2ORYB0AenimBVFZwqA5a2tBi+Xh6S+wxDl7HXVv1X7bX64FYV2FQ2PBFHabulETkEs0ye8PDw\nwBNPPIF33323xnp0ppw+fVqyXt7q1auxevVqTJo0Cc8//zy++eYbAKgR0qovZPzaa69BqVRiypQp\n+gWKN2/ezDXsHFBeWcWSJicypLNfg5up8NWIFmjjYZ0/E24qAWMC3TAm0A2iKCIptxzHbpUhMasM\n/8vR4FxO+V1b9Aj6AFQjEKkaFpycFZUByshrVVUtYk6Ku49Xqhhj076RvitERFJmmxUriqLJmai1\nGTRoENRqtdHrpq5VcnV1xZtvvok333yzTp9N9iWnVIdx32fi90xpV2Z3n4ru18ZsqTNFEAR09XFC\nVx8nABWLIutEETcLtbheqEVqQcXXm0Va5JTqkFOqQ3apDtklOuSW6VByZzB/Y6neheeiBJzvBBzn\nO0HISVHVmuR8JzSVFOTD39dbH5gqQ1JtXYQuSlQLZZXnUS1YCfoux8YYV0hE1NQ1ONjl5eXh2LFj\nePvtt9GzZ09z1ERUJ5klWjzyXRbOGoxPC23uhPjhzdHcwmPqGkohCGjrqUJbTxXCA+5+v06smGhQ\nUi6iWAuUlIvQQYRORNV/d+4T74yjU+iXBRCqLTMgQHlnbafK5Q6UggClUBnm6jebLjk5C8HBvnX/\nRhARUYPJDna+vr5G/5EXRRFt27bF2rVrzVYYkRy3irR45LtMXFCXS87383fCroda2OW2XgpBgLtK\ngHuTWYWSiIgai+wfDYsXL64R7ARBgI+PD+655x4MGzZMv38rUWO4XlCOMQcycTlf2jcZEeCMnQ81\nh5eT/YU6IiIiU2QnsZiYGEvWQVQnV/MrQt21AmmoG9rGBdsf9IO7iqGOiIgcD5vYqMm5lKvBmAOZ\nuFkkXTZkRDtXbBviB1d736+LiIjIiDoFu6tXr2L79u1ISUmBWq2GKEqXaBAEAV988YVZCySq7lyO\nBo98l4mMYmmoG9PBFVsG+8GZi2gREZEDkx3sdu7ciTlz5kCr1aJZs2bw9vaucY+97HVJtikhqwxj\nv8tCdqk01D3W0Q3vDvK16UVbiYiIGoPsYPfqq6+ic+fO+OSTTxAUFGTJmohq+L/bZRj3fSZyy6St\nxH8PdseGCB8oGeqIiIgge4R5Tk4Opk6dylBHje7YrVI8cqBmqHumqwc2DmCoIyIiqiQ72PXt2xep\nqamWrIWohh9vlmD8wSwUGGy5NfdeT6wJb8bdCIiIiKqRHexef/117N69G/Hx8Zash0jv4PUSPHYo\nC0UGoW5xLy+s6OvNMZ1EREQGZI+x69atG1544QU888wzmDt3Llq3bg2lUrpVkyAIOHHihNmLJMfz\nw40SPPlDFsqk8ySwrI83nuvpZZ2iiIiIbJzsYPfee+8hJiYGrq6u6NSpU62zYonM4cebtYe61+5v\nhtk9PK1TFBERURMgO9ht2LABYWFh2LFjB5o1a2bJmsiB/XKrFBMPZaNEuqEE1vVvhmldGeqIiIhM\nkT3GLj8/HxMmTGCoI4s5nl6Kxw9moVgrHVP3ZjhDHRERkRyyg92AAQNw9uxZS9ZCDuz322V47GAW\nCg0mSqy6vxme6cZQR0REJIfsYLdu3TqcOHEC69atQ0ZGhiVrIgdzUa3B+INZyNdIQ90rfb0xh2Pq\niIiIZJM9xq5v374QRRGrVq3CqlWr4OTkBIVCmgsFQcDNmzfNXiTZrxuFWjz6fc1twpbe5415IZz9\nSkREVBeyg93YsWO5bhiZVU6pDuO+z8T1QulMied6emJhKEMdERFRXckOdnFxcZasgxxMUbkOjx/M\nwgV1ueT85M7ueOk+LqVDRERUH7LG2BUVFWH06NH49NNPLV0POQCdKGLWzzk4dbtMcj6qvSve6u/D\nlmEiIqJ6khXs3N3dkZCQAK1We/ebie7itdP5+OpqieTcgFbO+HCwH1QKhjoiIqL6kj0rNiIiAseO\nHbNkLeQAdv5ZhLUJ+ZJz3XxU+OzB5nBVMdQRERE1hOxgt2bNGvz222946aWXcPXqVeh0uru/iKia\nk+mlmPtLjuRcC1cFdvy1OZo5y/6jSEREREbInjxx//33QxRFxMbGIjY2FgqFAk5OTpJ7uNwJGXOr\nSIvJR7Il+7+6KIHtw/zQwUv2H0MiIiIygcudkMVpdCKm/Dcb6cXSVt53BvgiLMDFSlURERHZHy53\nQha34rc8HE+XzoD9V09PTOjkbqWKiIiI7BMHNpFF7b1ajE1nCyTnhrRxwQu9uVYdERGRudUp2OXk\n5GDZsmUIDw9HmzZt0KZNG4SHh2P58uXIycm5+xuQQ0ktKMc/jkr/XPzFXYktg32h5LImREREZic7\n2F2/fh2DBg3Cpk2b4ObmhtGjR2P06NFwd3fHxo0bMWjQIFy/ft2StVITotWJePbnHOSVifpzTgrg\n46F+aOGqtGJlRERE9kv2GLvly5cjNzcX+/btw4ABAyTXjh07hokTJ+KVV17BBx98YPYiqel5538F\nOHpLOq7ulb7N0K+ls5UqIiIisn+yW+wOHz6MmTNn1gh1QMXixTNmzMAPP/xg1uKoabpQIGDl73mS\nc8PauODZ7h5WqoiIiMgxyA52xcXFaNGihdHrLVq0QHFxsVmKoqarTCvi5Ysu0FRb2cTXRUDsIF8o\nuFwOERGRRckOdl27dsWuXbtQWlpa41pZWRm++OILdOvWzazFUdOz/kw+LhdJ/1htjPBFa3eOqyMi\nIrI02WPs5s+fj6effhpDhw7F1KlTERQUBABITk7GRx99hAsXLmDbtm0WK5Rs3wW1psY+sJOC3DEm\n0M1KFRERETkW2cEuOjoamzdvxrJly7Bo0SL9LhSiKKJly5aIi4vD6NGjLVYo2TadKGLeUbWkC9bf\nVYHX7m9mvaKIiIgcjNFgl5eXBw8PDyiVVV1ojz/+OMaNG4fTp08jNTUVoiiiffv26N27N1Qq7vfp\nyD5OKsLJDOks2DXhzeDrwjWwiYiIGovRn7qBgYGIj4/XH8+ZMwf/93//B5VKhX79+uHRRx/FuHHj\n0K9fv3qHuqNHj2LixIno1q0bfHx8sH37dsl1URSxevVqdO3aFa1atUJUVBTOnz8vuUetVmPGjBlo\n37492rdvjxkzZkCtVterHqqfnFIdXv09V3JuZDtXPMIuWCIiokZlNNg5OztLJkp89tlnuHLlilk/\nvLCwEN27d8frr78ON7eaIWDjxo2IjY3FG2+8gcOHD8Pf3x9jx45Ffn7VOK7p06cjMTERu3btwu7d\nu5GYmIiZM2eatU4y7bXTecgprVqI2E0hYm14M313PRERETUOo01tXbp0wdatW9GiRQt4eXkBAJKS\nknD06FGTb1jbOnfGDB8+HMOHDwcAzJ49W3JNFEXExcVh/vz5iI6OBgDExcUhODgYu3fvxpQpU5CU\nlIRDhw7hwIEDCAsLAwCsX78ekZGRSE5ORnBwsOxaqH7O5Wiw9UKh5NzUdhq09WTXPBERUWMz+tN3\nxYoVePrppzFp0iQIggBBEPDWW2/hrbfeqvV+URQhCAKys7PNUlhKSgrS09MxbNgw/Tk3NzdERETg\n5MmTmDJlCk6dOgVPT099qAOA8PBweHh44OTJkwx2FiaKIp4/mQttVWMdAr2UmPSXcusVRURE5MCM\nBrvBgwfjzJkzOHv2LDIyMvDUU09h5syZ6N+/f6MUlp6eDgDw9/eXnPf390daWhoAICMjA82bN5d0\n+QmCgBYtWiAjI8PoeycnJ1ug4sb/DGs7mq3AT2muknP/aFsEF4VjPL8pfH4+vyPj8/P5HZmln/9u\njVYm+8s8PT0RHh4OoKKLNTIyEoMHDzZfdTIYjtOqbBk0dr22ewxZuiXPEbqBdaKIqXtvA9Dozw1p\n44JpYW1w6dIlu39+Uxzh998UPj+fn8/P53dUtvD8sgdC7du3z5J11BAQEACgolWubdu2+vOZmZn6\nVryWLVsiMzNTEuREUURWVlaNlj4yr6+uFuNMtkZy7pW+3pwwQUREZEX1XmQsIyMDfn5++PHHH81Z\nj16HDh0QEBCAI0eO6M+VlJTg+PHj+jF1999/PwoKCnDq1Cn9PadOnUJhYaFk3B2ZV7lOxKrfpTtM\nPHqPG0KbO1upIiIiIgLq0GJXG1EU736TCQUFBbh8+TIAQKfT4fr160hMTISvry/atWuHWbNmYd26\ndQgODkZQUBDWrl0LDw8PjB8/HkDFzN2//vWvWLBgATZu3AhRFLFgwQKMGDHC6k2h9uyzS0W4lFc1\nQUIpAC/09rJiRURERAQ0MNg11OnTpyXbkK1evRqrV6/GpEmTEBcXh3nz5qG4uBiLFi2CWq1Gnz59\nEB8fr19+BQA++OADLFmyBI8++igAIDIyEmvWrGn0Z3EU5ToR6xOlrXVPBLkjqJmTlSoiIiKiSlYN\ndoMGDTK5S4QgCIiJiUFMTIzRe3x9ffH+++9bojyqxVdXi3ElX6s/dlIAi3uxtY6IiMgW1HuMnbe3\nN2JjY9GtWzdz1kM2TBRFrD9TIDn3eCd3tONixERERDah3j+RXV1d8cQTT5izFrJxB6+X4my1mbAC\ngHkhntYriIiIiCSMBrvU1NR6vWG7du3qXQzZtvVnpGPrxgS6Iphj64iIiGyG0WDXs2fPeq1JZq4t\nxci2/JFZhuPpZZJzC0I4to6IiMiWGA1277zzDhebJb0tFwolx4Nbu6BXC65bR0REZEuMBrsnn3yy\nMesgG5ZdosXuy0WSczO6eVipGiIiIjKmXpMnKrftAoDmzZuzZc/O/Tu5CCVVK5ygnacSI9u5Wq8g\nIiIiqlWdlju5fPkynn76abRv3x6dO3dG586d0b59e0ydOlW/gwTZF61OrNENO72rB5QKhnkiIiJb\nI7vF7vz58xgxYgRKSkoQGRmJzp07AwAuXryI/fv34/Dhw/j222+5rp2dOXSjFKkFVc11rkrg78Hu\nVqyIiIiIjJEd7JYvXw53d3f897//RceOHSXXrly5gsjISLzyyivYsWOH2Ysk69meLG2tG9fRHX6u\nSitVQ0RERKbI7oo9fvw4pk+fXiPUAcA999yDadOm4dixY2Ytjqwrq0SLb1NLJOcms7WOiIjIZskO\ndlqtFi4uLkavu7q6QqvVGr1OTc+uy8XQ6KqOg5upcH9LLnFCRERkq2QHu9DQUHzyySdQq9U1rqnV\nanzyySfo1auXWYsj69qeLF3i5Ikgd86AJiIismGyx9i98MILGDt2LPr27YsnnngCQUFBAIDk5GTs\n2LEDubm52LBhg8UKpcaVmFWGM9X2hVUIwMQgdsMSERHZMtnBbuDAgfjPf/6DpUuXYtOmTZJroaGh\n2Lp1KwYMGGD2Ask6vvizWHL8YBsXtHbnpAkiIiJbVqcFih944AH89NNPSE9PR2pqKkRRRPv27REQ\nEGCp+sgKdKKIL69Kgx1b64iIiGxfvXaeCAgIYJizY/93uwzXC6smwrirBO40QURE1ATUKdiJooif\nfvoJKSkpUKvVEEVRcl0QBPzzn/80a4HU+OKvSFvrRrR1hYdTnTYpISIiIiuQHez++OMPTJkyBSkp\nKTUCXSUHU9KSAAAgAElEQVQGu6ZPJ4r4yqAb9pF73KxUDREREdWF7GA3b948ZGdnY/369ejTpw+8\nvb0tWRdZyYn0MqQVVS1e56kSMLwtu2GJiIiaAtnBLikpCS+++CKeeuopS9ZDVrY3RdpaF9neFW4q\nrl1HRETUFMgeONWpUyejXbBkH0RRxLfXpFuIjQlkNywREVFTITvYxcTE4IMPPkBqaqol6yEruqAu\nR0pB1WxYFyUwrI3xbeSIiIjItsjuih01ahRKSkoQFhaGgQMHok2bNlAqpQvWCoKAtWvXmr1IahwH\nUqWtdQ+0cuFsWCIioiZEdrD76aefsGDBAhQXF+PgwYO13sNg17QZBruR7TlpgoiIqCmRHeyef/55\neHt745NPPuGsWDt0u1iLUxllknMjOBuWiIioSZHdz3blyhX885//xNChQxnq7ND310tQfWpMiJ8T\n2nrWa2MSIiIishLZwa5r167IycmxZC1kRYeul0qOuYUYERFR0yM72L366qvYtm0bTp48acl6yAp0\noogf06TBjosSExERNT2y+9o2btwIDw8PREZGIigoCG3btq11VuwXX3xh9iLJshKzNMgurdptwttZ\nQO8WTlasiIiIiOpDdrC7cOECBEFA27ZtUVJSgkuXLtW4RxC4Q0FTZNha90ArF6gU/L0kIiJqamQH\nuzNnzliyDrKiIzelwW4IFyUmIiJqkkyOsVu7di0DnZ0rLhdxPF0a7Ia24fg6IiKipshksNu6dSsG\nDx6M7t27Y968efjmm29QWFjYWLVRIziZUYrSql3E0NZDiY7eSuMvICIiIptlsiv23LlzSEhIwPff\nf4+DBw/i008/hZOTEyIiIjBixAgMHz4c99xzT2PVShbwU5pha50Lx0oSERE1UXdd7iQ0NBSLFi3C\n999/j+TkZGzYsAG+vr5YvXo1+vTpg379+uHFF1/ETz/9hPLy8saomczo2C3pbhODWnN8HRERUVNV\npx3e/fz8MHHiRHz44Yf4888/8fXXXyMqKgpHjhxBdHQ0OnbsiMmTJ+PUqVOWqpfMqKRcxO+Z0mAX\nEeBspWqIiIiooeoU7KpTKpUYMGAAli9fjmPHjiExMRHLli1DaWmp2RYx1mq1WLlyJXr27ImAgAD0\n7NkTK1eulLQMiqKI1atXo2vXrmjVqhWioqJw/vx5s3y+vfs9swxlVcvXoZ2nktuIERERNWFm+yne\nrl07TJ8+HdOnTzfXW2LDhg3YsmUL4uLi0L17d/zvf//DrFmz4OzsjMWLFwOoWDg5NjYWsbGxCA4O\nxpo1azB27Fj8+uuv8PLyMlst9uh4urS1rj9b64iIiJo0k8Fuz549dXozhUIBb29vdOnSBW3atGlQ\nYQBw6tQpjBw5EpGRkQCADh06IDIyEr/99huAita6uLg4zJ8/H9HR0QCAuLg4BAcHY/fu3ZgyZUqD\na7BnhsucRARwfB0REVFTZjLYTZ06FYIgQBTFOr2pIAh44IEHsGXLFrRo0aLexYWHh+PDDz/ExYsX\n0blzZ1y4cAE///wzFixYAABISUlBeno6hg0bpn+Nm5sbIiIicPLkSQY7E7Q6Eacy2GJHRERkT0wG\nu6+//rpObyaKIgoKCvDbb7/hnXfeweLFi7F169Z6Fzd//nwUFBQgLCwMSqUS5eXlWLhwob67Nz09\nHQDg7+8veZ2/vz/S0tKMvm9ycnK9a5KrMT6jIZIKBORp3PTHzVQihIyrSL5tnve39ee3ND4/n9+R\n8fn5/I7M0s8fHBxs8rrJYDdw4MB6fWhkZCREUcSWLVvq9fpK8fHx2LFjB7Zs2YKuXbvizJkzeP75\n59G+fXtMnjxZf5/humuiKJpci+1u35SGSk5OtvhnNNSRcwUAcvXHA9q4oXPntmZ576bw/JbE5+fz\n8/n5/I6Kz2/957fYFMiBAwc2eNmTZcuW4R//+AfGjRsHAOjRowdSU1Oxfv16TJ48GQEBAQCAjIwM\ntG1bFUoyMzNrtOKRVEK2RnIc1pLdsERERE2d0eVOVq9eDbVaXec3zMnJwWuvvYahQ4fWuSvXUFFR\nEZRK6fZWSqUSOl3FGh0dOnRAQEAAjhw5or9eUlKC48ePIywsrEGfbe8SsqTBrldzJytVQkREROZi\nNNjt27cP9957L+bMmYNDhw6htLTU2K0oKSnBoUOHMHv2bISEhODbb781S3EjR47Ehg0b8N133yEl\nJQVff/01YmNjMWrUKAAVXbCzZs3Chg0bsHfvXpw7dw6zZ8+Gh4cHxo8fb5Ya7FFJuYgLOdJg17M5\nW+yIiIiaOqNdsUePHsWuXbuwadMmfPbZZ1CpVOjSpQsCAwPh4+MDURShVquRkpKCpKQkaLVa9OzZ\nExs3btR3nTbUmjVrsGrVKvzrX/9CZmYmAgIC8NRTT+nXsAOAefPmobi4GIsWLYJarUafPn0QHx/P\nNexMOK/WoLzaROd2nkr4utR7rWoiIiKyESbH2E2YMAETJkxAQkIC9u/fj19//RV//PEHsrOzAVRs\nMdalSxeMGTMGDz/8MO69916zFufl5YXXX38dr7/+utF7BEFATEwMYmJizPrZ9izRoBs21I/dsERE\nRPZA1uSJ0NBQhIaGWroWaiSG4+tCOb6OiIjILrD/zQElZEkXJg7l+DoiIiK7wGDnYMp1Iv5XY+IE\nW+yIiIjsAYOdg7mYW44SbdVxgJsCrdyVxl9ARERETQaDnYPh+DoiIiL7xWDnYAzH13H9OiIiIvvB\nYOdgDFvsenKpEyIiIrshO9itWLECFy9etGQtZGE6UcTZbHbFEhER2SvZwe6dd95BeHg4hgwZgs2b\nN+P27duWrIss4EqeFvmaqi0nfJwFtPfkxAkiIiJ7ITvYJSUlYc2aNXBxcUFMTAy6d++Oxx57DPHx\n8SgpKbFkjWQmta1fJwiClaohIiIic5Md7Hx9fTF9+nR89913OH36NBYuXIgrV65g2rRp6Ny5M+bO\nnYuff/7ZkrVSA3FGLBERkX2r1+SJwMBALFmyBL/++it++OEHDB06FNu3b0d0dDTuvfderFy5Emlp\naeaulRooMZsLExMREdmzes+KLSwsxI4dO/Dqq69i3759UKlUGDlyJPr27YtNmzahT58++Oqrr8xZ\nKzWAKIpssSMiIrJzqrrcrNPpcPjwYezcuRPffPMNioqKEBoaitdeew0TJkyAn58fAOD27duYOnUq\nli5diujoaIsUTnVzvVCL7FKd/thTJaCTd51++4mIiMjGyf7JHhMTg/j4eNy+fRsBAQGYPn06Jk2a\nhK5du9a419/fH3/729/w7LPPmrVYqj/D1rqQ5k5QcOIEERGRXZEd7LZt24aoqChMnDgRQ4cOhUJh\nuhc3PDwcsbGxDS6QzKNGsOPCxERERHZHdrBLSkqCl5eX7Dfu0KEDOnToUK+iyPwMJ05wfB0REZH9\nkT15YuDAgfjmm2+MXj9w4ABCQ0PNUhSZX2Ita9gRERGRfZEd7K5du4bCwkKj1wsLC5GammqWosi8\n0ou0SCuqmjjhogS6+HDiBBERkb2p03InpnYpuHTpUp26aqnxGHbDdvd1gpOCEyeIiIjsjclmm88+\n+wyff/65/njt2rXYtm1bjfvUajXOnTuHkSNHmr9CarBEw/XrOHGCiIjILpkMdsXFxcjKytIfFxQU\n1Dob1sPDA1OnTsWSJUvMXyE1WG17xBIREZH9MRnspk2bhmnTpgEAevbsiddffx0PP/xwoxRG5sMd\nJ4iIiByD7BH0iYmJlqyDLERdqkNKgVZ/rBQqxtgRERGR/an3XrHUNBhOnOjio4KrihMniIiI7JHR\nFjtfX18oFAqkpaXB2dkZvr6+JmfFAhWzZquPySPr4/g6IiIix2E02C1evBiCIEClUkmOqWmpMSOW\n4+uIiIjsltFgFxMTY/KYmgZOnCAiInIcHGNnxwo1OiTnluuPBQD3cg07IiIiuyU72K1fvx4jRoww\nej0yMhKbNm0yS1FkHmezNRCrHXfyVsHLiVmeiIjIXsn+Kb9r1y7069fP6PV+/fphx44dZimKzIPd\nsERERI5FdrC7evUqgoODjV7v1KkTUlJSzFIUmUeCwVInPRnsiIiI7JrsYOfi4oJbt24ZvZ6Wllbr\ndmNkPZwRS0RE5FhkJ7H7778fn376KXJzc2tcU6vV+PTTTxEWFmbW4qj+SrUiLqgNWuw4cYKIiMiu\nyd5S7Pnnn0dkZCQGDBiAWbNmoVu3bgCAc+fOYfPmzbh9+zY+/vhjS9VJdXQ+RwONruq4rYcSfq5K\n6xVEREREFic72PXu3Rs7d+7EvHnzsHTpUv1ixaIoIjAwEDt37kTfvn0tVijVjeFWYhxfR0REZP9k\nBzsAGDx4ME6fPo2EhARcuXIFoiiiY8eOCA0N5a4UNsZwfB27YYmIiOxfnWc7CIKAXr16YezYsXj0\n0UfRq1cvi4a6W7du4dlnn0WnTp0QEBCAsLAw/PLLL/rroihi9erV6Nq1K1q1aoWoqCicP3/eYvU0\nFZw4QURE5Hjq1GIHABqNBhcvXkReXh50Ol2N6wMGDDBLYUDFpIwRI0YgPDwcX3zxBZo3b46UlBT4\n+/vr79m4cSNiY2MRGxuL4OBgrFmzBmPHjsWvv/4KLy8vs9XSlGh1Is7mGHbFOlupGiIiImossoOd\nKIp49dVX8cEHH6CwsNDofdnZ2WYpDADefvtttGrVCu+9957+XGBgoKSmuLg4zJ8/H9HR0QCAuLg4\nBAcHY/fu3ZgyZYrZamlKLuWVo6i8as+JFq4KtHHnUjRERET2TvZP+w0bNmD9+vUYN24cNm/eDFEU\nsXz5cqxfvx7dunVDSEgI9uzZY9bi9u/fjz59+mDKlCkICgrCwIED8f7770MUK0JLSkoK0tPTMWzY\nMP1r3NzcEBERgZMnT5q1lqaktvF1HANJRERk/2S32H366acYM2YMNmzYoG+VCw0NxeDBgzFx4kQ8\n+OCD+OWXXzB48GCzFXf16lV8+OGHmD17NubPn48zZ85gyZIlAIAZM2YgPT0dACRds5XHaWlpRt83\nOTnZbDVa8zOM+e8VJwBVY+raKQqQnJzTqDVY8/ltAZ+fz+/I+Px8fkdm6ec3tQsYUIdgd/36dcyZ\nMwcA9DtMlJaWAqjYleLxxx/He++9hxdffLG+tdag0+nQu3dvvPzyywAqguTly5exZcsWzJgxQ3+f\nYWuUKIomW6ju9k1pqOTkZIt/himpf2YCKNUfDwlqieB73Bvt8639/NbG5+fz8/n5/I6Kz2/955fd\nFevj44OSkhIAgLe3N5ydnXHjxg39dRcXF7OOrwOAgIAAdOnSRXKuc+fOuH79uv46AGRkZEjuyczM\nrNGK5yhEUURCVpnkXE8/TpwgIiJyBLKDXbdu3XDmzJmKFykUuO+++/Dhhx/ixo0bSE1Nxccff2z2\nlBoeHo5Lly5Jzl26dAnt2rUDAHTo0AEBAQE4cuSI/npJSQmOHz/usNubXSvQIresauKEl5OAe7y5\n4wQREZEjkB3sJkyYgKSkJH2r3bJly3Dp0iWEhIQgNDQUf/75J5YtW2bW4mbPno1ff/0Va9euxeXL\nl/Hll1/i/fffx/Tp0wFUdMHOmjULGzZswN69e3Hu3DnMnj0bHh4eGD9+vFlraSoMd5y4188JCk6c\nICIicgiyx9g9+eSTePLJJ/XH/fv3x4kTJ/DNN99ApVLhwQcfRKdOncxa3H333Yft27djxYoVePPN\nN9G2bVu88MIL+mAHAPPmzUNxcTEWLVoEtVqNPn36ID4+3mHXsEvgwsREREQOq84LFFcXGBiI2bNn\nm6uWWo0YMQIjRowwel0QBMTExCAmJsaidTQVZ2qMr2OwIyIichR1DnaXL1/G999/j9TUVABA+/bt\n8dBDD6Fjx45mL47qzrDFjjtOEBEROQ7ZwU6r1WLJkiX46KOPamwl9sILL+Dpp5/GmjVroFRyoL61\npBdpcau46vfGRQl08WlQoywRERE1IbInT6xcuRIffvghHnvsMRw5cgTXrl3DtWvXcOTIEUyYMAEf\nffQRVq5cacla6S4MJ05093WCk4ITJ4iIiByF7Oac7du345FHHkFcXJzkfK9evbB582YUFxdj+/bt\n+sWEqfEZbiUWyvF1REREDkV2i11RUREGDhxo9PoDDzyAoqIisxRF9ZOYbTBxguPriIiIHIrsYBcR\nEYETJ04YvX7ixAlERESYpSiqHy51QkRE5NhkB7t169YhMTER//rXv5CUlASNRgONRoOkpCQ899xz\nOHPmDN566y1L1komqEt1uJqv1R8rhYoxdkREROQ4ZI+x69evH0RRxMWLF/HRRx9BuLObgShWbF+l\nUqnQr18/yWsEQcDNmzfNWC4Zc8Zg4kTnZiq4qThxgoiIyJHIDnZjx47VhzmyPYYzYnuyG5aIiMjh\nyA52hrNhybYkGu44wYkTREREDkf2GDuybTWWOmGLHRERkcOp87YEN2/eREJCAvLy8mrsQAEAkyZN\nMkthJF9RuQ5JueWScyFcw46IiMjhyA52paWlmDNnDvbs2QOdTgdBEPQTJ6qPvWOwa3zncsqhE6uO\nA72UaObMxlgiIiJHI/un/6pVq/Dll1/ixRdfxL59+yCKIuLi4rBnzx4MGzYMISEhOHr0qCVrJSPY\nDUtERERAHYLdnj17MHHiRDz33HPo1q0bAKB169YYMmQIdu3aBXd3d2zdutVihZJxCYYTJ/w4cYKI\niMgRyQ52GRkZ+nXqVKqKHtySkhIAFV2x0dHR2Lt3rwVKpLsxXOqELXZERESOSXawa9GiBfLy8gAA\nXl5ecHNzw5UrV/TXNRoNCgsLzV8hmaTRiTiXwzXsiIiIqA6TJ0JCQvDbb78BqGihGzBgAOLi4hAa\nGgqdTof3338fISEhFiuUapekLkdp1U5iaO2uQEs3pfUKIiIiIquR3WL31FNPoby8XN/9umLFChQU\nFCAqKgqjRo1CUVERVq1aZbFCqXY1FibmMidEREQOS3aLXWRkJCIjI/XH3bp1w+nTp/Hzzz9DqVQi\nPDwcPj4+FimSjEswmBEbwh0niIiIHFadFyiuztvbG1FRUeaqheqBEyeIiIioksmu2NzcXIwbNw7r\n1q0z+SZr167F+PHjUVBQYNbiyDSdKOKsQbBjVywREZHjMhnstmzZglOnTmHy5Mkm32Ty5Mk4efIk\n17FrZFfytMjXVG054eMsoL0nJ04QERE5KpPBbv/+/XjkkUfg7+9v8k1atmyJRx99FF999ZVZiyPT\nErMNJk40d5Zs70ZERESOxWSwu3jxInr37i3rjUJDQ3Hx4kWzFEXyGE6cYDcsERGRYzMZ7LRarX6X\nibtRqVQoLy83S1EkD/eIJSIioupMBrs2bdrg7Nmzst7o7NmzaN26tVmKorsTRbFmix2DHRERkUMz\nGeyGDRuGnTt3Ij093eSb3Lp1Czt37sSDDz5o1uLIuJtFOmSV6vTH7ioBQd4NWr2GiIiImjiTwe6f\n//wntFotxowZg1OnTtV6z6lTpxAdHQ2tVou5c+dapEiqyXDHiXt9naBUcOIEERGRIzPZxNOuXTts\n27YNU6ZMwciRI9GhQwf06NEDnp6eKCgowLlz53D16lV4eHjgo48+Qvv27RurbofHblgiIiIydNe+\nuwcffBC//PILNm7ciAMHDmD//v36a61bt8aUKVMwd+5cBAYGWrJOMmC44wSDHREREckalNW+fXus\nW7cO69atQ35+PvLz8+Hl5QUvLy9L10dGGM6I5VInREREVOfR9gx01pddosX1Qq3+2EkBdPNlsCMi\nInJ0JidPkG0y7Ibt6uMEFyUnThARETk6BrsmyHDiBBcmJiIiIoDBrkni+DoiIiKqDYNdE8SlToiI\niKg2soNdamoqiouLjV4vLi5GamqqWYoyZt26dfDx8cGiRYv050RRxOrVq9G1a1e0atUKUVFROH/+\nvEXrsKZ8jQ5/5lXtySsAuJctdkRERIQ6BLvQ0FDs27fP6PVvv/0WoaGhZimqNr/++iu2bduGHj16\nSM5v3LgRsbGxeOONN3D48GH4+/tj7NixyM/Pt1gt1nQ2WwOx2nFQMxU8ndjwSkRERHUIdqIomrxe\nXl4OQbDMzMzc3Fw888wz2LRpE3x8fCQ1xcXFYf78+YiOjkb37t0RFxeHgoIC7N692yK1WJvh+DpO\nnCAiIqJKdWrqMRbccnNzcejQIfj7+5ulKEOVwW3w4MGS8ykpKUhPT8ewYcP059zc3BAREYGTJ09a\npBZrq7HjBLthiYiI6A6TCxS//vrrWLNmDYCKUDdjxgzMmDGj1ntFUcTs2bPNXuC2bdtw+fJlvPfe\nezWupaenA0CNQOnv74+0tDSj75mcnGzeIhvxM07ddEX1PO5XkoHk5FsW+ayGaIzvsS3j8/P5HRmf\nn8/vyCz9/MHBwSavmwx2ffr0wbRp0wAAW7ZswdChQ9GpUyfJPYIgwMPDA7169cKYMWMaWK5UcnIy\nVqxYgW+//RbOzs5G7zNsSRRF0WS38N2+KQ2VnJxskc8o1Yq4cvSm5NzDIYHwc1Wa/bMawlLP31Tw\n+fn8fH4+v6Pi81v/+U0Gu4ceeggPPfQQAKCwsBBTp05F3759G6UwADh16hSysrLQv39//TmtVotj\nx45h69atOHHiBAAgIyMDbdu21d+TmZlpsW5hazqfo0F5taGObT2UNhfqiIiIyHpk7xX77rvvWrKO\nWkVFRaF3796Sc3PmzEGnTp3w3HPPISgoCAEBAThy5Ajuu+8+AEBJSQmOHz+OFStWNHq9lmY4vo4T\nJ4iIiKg62cGu0s2bN5GQkIC8vDzodLoa1ydNmmSWwgDAx8dHMgsWANzd3eHr64vu3bsDAGbNmoV1\n69YhODgYQUFBWLt2LTw8PDB+/Hiz1WErauw4wWBHRERE1cgOdqWlpZgzZw727NkDnU4HQRD0S6BU\nH89mzmAnx7x581BcXIxFixZBrVajT58+iI+Ph5eXV6PW0RgSssokx2yxIyIioupkB7tVq1bhyy+/\nxIsvvoiwsDCMGjUKcXFxaNWqFd555x3cvn0bmzdvtmStAID9+/dLjgVBQExMDGJiYiz+2dak1Yk4\nm10uOdfTz/iEEiIiInI8stex27NnDyZOnIjnnnsO3bp1AwC0bt0aQ4YMwa5du+Du7o6tW7darFBH\nl5xXjmJt1cyJFq4KtHbnjhNERERURXYyyMjIQL9+/QAAKlVFQ19JSQmAilaz6Oho7N271wIlElD7\njhOW2umDiIiImibZwa5FixbIy8sDAHh5ecHNzQ1XrlzRX9doNCgsLDR/hQQASDCcOMEdJ4iIiMiA\n7DF2ISEh+O233wBUtNANGDAAcXFxCA0NhU6nw/vvv4+QkBCLFeroEmtMnOD4OiIiIpKS3WL31FNP\noby8XN/9umLFChQUFCAqKgqjRo1CUVERVq1aZbFCHZkoijX3iOWMWCIiIjIgu8UuMjISkZGR+uNu\n3brh9OnT+Pnnn6FUKhEeHl5jzTkyj5QCLXLLqiZOeDkJCPTijhNEREQkVecFiqvz9vZGVFSUuWoh\nIwwnToT4OUHBiRNERERkoE7rZZSVleGTTz7BM888g0ceeQQJCQkAALVajc8//xw3btywSJGOjjtO\nEBERkRyyW+yys7MxevRonDt3Di1btsTt27ehVqsBVLTcrVq1ChcuXMArr7xisWIdVWI2J04QERHR\n3clusXv55ZeRmpqKAwcO4NixY/rtxABAoVBgzJgxOHjwoEWKdHQ1Wuy41AkRERHVQnawO3DgAGbO\nnImwsLBaF8bt1KkTrl+/btbiCEgv0uJWsU5/7KoEuvg0aGgkERER2SnZwS4/Px9t27Y1er20tBRa\nrdYsRVEVw2VOuvs6QaXgxAkiIiKqSXaw69ixI06fPm30+uHDh/V7yJL5cMcJIiIikqtOCxR/9tln\n+OKLL6DTVXQNCoKAoqIiLF++HIcPH8aUKVMsVqij4o4TREREJJfswVozZ87EhQsXMHPmTHh5eQEA\npk6dCrVaDa1Wi+nTp+PJJ5+0WKGOqkaLHZc6ISIiIiPqNAp//fr1mDhxIvbs2YPLly9Dp9Phnnvu\nwdixYxEREWGpGh2WulSHlIKqcYtKoWKMHREREVFt6jy9MiwsDGFhYZaohQycMZg40aWZCm4qTpwg\nIiKi2tVp5wlqXAkG4+tC2A1LREREJphssRs9enSd3kwQBOzdu7dBBVEVw6VOOHGCiIiITDEZ7DIz\nMyWLEWu1WiQnJ6Ndu3bw8PCweHGOjnvEEhERUV2YDHbHjx+XHGdlZSEoKAhvv/02Bg8ebNHCHF1R\nuQ4Xc8sl50K4hh0RERGZUKcxdrVtJUaWcS6nHLqq7Xhxj5cSzZw5JJKIiIiMY1KwUYYTJzi+joiI\niO6Gwc5GcXwdERER1RWDnY3iHrFERERUVyYnT9y+fVtynJ2dDQDIzc2tca2Sv7+/mUpzXBqdiHM5\nbLEjIiKiujEZ7Dp37lzrhImnn37a6Gsqwx/V3wV1Ocp0Vcet3RVo6aa0XkFERETUJJgMdosXL+ZM\nWCtINJg40ZMTJ4iIiEgGk8EuJiamseqgampMnOD4OiIiIpKBkydskOFWYhxfR0RERHIw2NkYnSji\nTJbhHrEMdkRERHR3DHY25nJeOQrKq7ac8HEW0M6DEyeIiIjo7hjsbIzh+LrQ5s6cwEJERESyMNjZ\nGI6vIyIiovpisLMx3HGCiIiI6kt2sDtx4gTWr19v9Pr69etx6tQpsxTlqERRrKUrlsGOiIiI5DG5\njl11b7zxBnx8fIxeP3v2LH755Rf85z//MUthjuhGoRZZpVVbTrirBHTylv1bRERERA5OdotdYmIi\n7r//fqPX+/Xrh4SEBLMUVemtt97C0KFD0a5dO3Tq1AmPP/44zp07J7lHFEWsXr0aXbt2RatWrRAV\nFYXz58+btY7GYji+LsTPCUoFJ04QERGRPLKDXVFR0V1nZxYUFDS4oOp++eUXTJs2Dd999x327t0L\nlXs+K7oAACAASURBVEqFRx55BDk5Ofp7Nm7ciNjYWLzxxhs4fPgw/P39MXbsWOTn55u1lsbA8XVE\nRETUELKDXVBQEA4fPmz0+qFDh9CxY0ezFFUpPj4ef/vb39C9e3f06NED7733HjIzM3HixAkAFa11\ncXFxmD9/PqKjo9G9e3fExcWhoKAAu3fvNmstjaHGVmIcX0dERER1IHsA1+TJk7F48WIsXrwYMTEx\n8PX1BQBkZ2dj9erVOHz4MFatWmWxQoGKFkGdTqcf65eSkoL09HQMGzZMf4+bmxsiIiJw8uRJTJky\npdb3SU5Otmid9f2M3zNcUT1r+xamITlZNP4CG9YY32Nbxufn8zsyPj+f35FZ+vmDg4NNXpcd7J55\n5hmcOXMGH3zwAbZs2YKWLVsCADIyMiCKIp544gnMmjWrYdXexfPPP4+QkBD9WL/09HQAgL+/v+Q+\nf39/pKWlGX2fu31TGio5ObnOn5FVokX6L7f0x04KYHjPTnBWNr0xdvV5fnvC5+fz8/n5/I6Kz2/9\n56/TlMu3334bEyZMwN69e3H16lWIooh77rkH0dHRGDhwoKVqBAC88MILOHHiBA4cOAClUrrFluHY\nP1EUm9xuDYbdsN18nJpkqCMiIiLrqfNaGoMGDcKgQYMsUYtRMTExiI+Px9dff43AwED9+YCAAAAV\nrYZt27bVn8/MzKzRimfrakyc4Pg6IiIiqiOb33liyZIl2L17N/bu3YvOnTtLrnXo0AEBAQE4cuSI\n/lxJSQmOHz+OsLCwxi61QQyXOuHCxERERFRXRlvsRo0aBYVCgfj4eKhUKowePfqubyYIAvbu3Wu2\n4hYuXIidO3fi008/hY+Pj35MnYeHBzw9PSEIAmbNmoV169YhODgYQUFBWLt2LTw8PDB+/Hiz1dEY\nasyI5VInREREVEdGg50oitDpqnZB0Ol0dx23JormncG5ZcsWAEB0dLTk/JIlSxATEwMAmDdvHoqL\ni7Fo0SKo1Wr06dMH8fHx8PLyMmstlpSv0eFSXrn+WADQg8GOiIiI6shosNu/f7/J48agVqvveo8g\nCIiJidEHvaborEE3bHAzFTydbL6XnIiIiGyMrPRQXFys39mBzI8TJ4iIiMgcZAU7Nzc3rF+/Htev\nX7d0PQ7JcHxdKLthiYiIqB5k9/fde++9uHz5siVrcVgJWWWSY7bYERERUX3IDnYvvfQStm3bhu++\n+86S9TicUq2IJHW55FzP5s5WqoaIiIiaMtkLFL/zzjvw9fXFpEmT0KZNGwQGBsLNzU1yjyAI+OKL\nL8xepD07n6NBebXJxO08lfB14cQJIiIiqjvZwe7ChQsQBEG/w8O1a9dq3NPUtvGyBTUmTnB8HRER\nEdWT7GB35swZS9bhsLjjBBEREZkL+/ysjBMniIiIyFxkt9hVl5+fj7y8PMnOFJXatWvX4KIchVYn\n4n/Z0okToZw4QURERPVUp2D38ccfY9OmTbhy5YrRe7KzsxtclKNIzitHsbZq5oS/qwKt3NiISkRE\nRPUjO0V88sknWLBgATp06IClS5dCFEXMmjULCxYsQMuWLRESEoJNmzZZsla7U9uOE5yAQkRERPUl\nO9ht3rwZQ4YMQXx8PJ5++mkAwPDhw/HSSy/hxIkTUKvVyMvLs1SddqnGjhMcX0dEREQNIDvYXb58\nGQ8//HDFixQVL9NoKoKJj48PJk+ejC1btligRPtVY+KEH8fXERERUf3JDnYeHh4QxYrxYJ6enlAq\nlUhLS9Nf9/Pzw82bN81foZ0SRZFLnRAREZFZyQ52nTt3RlJSEgBApVIhJCQEO3fuhEajQUlJCXbu\n3IkOHTpYrFB7k1KgRV5Z1cQJbycBHbyUVqyIiIiImjrZwe7hhx/GgQMHUFJSAgBYuHAhjh07hsDA\nQAQFBeHkyZNYsGCBxQq1N4YTJ0KaO0HBiRNERETUALKXO5k7dy7mzp2rP46KisL+/fvx1VdfQaVS\nYeTIkRg4cKBFirRHZ7iVGBEREZlZvRYortS/f3/079/fXLU4lJo7TnDiBBER/X979x7X4/0/fvxR\nUUnyppOKhHKK5DDlOOFja4dQksN8KMk228yEGnMY+4TMufnwwfjsM2bLKcOwKcda8yEMmcNU87NM\nVhTR4f39w8/74+rd4Y13xz3vt5vbzfW6Xtf1fr5eV+/rer5f10mI5/PUid3du3c5cuQIaWlpADg6\nOtKzZ0/q1aun9+BqKrVaTZI86kQIIYQQevZUid3y5ctZuHAh9+7d09whC2BmZsbUqVOZOHGi3gOs\niW7cK+SP3P+9jq2OkQEt6z/X4KkQQgghhO6J3YoVK5g1axY9e/YkODgYZ2dn1Go1V65cYe3atcyZ\nMwdDQ0PFdXiieElFTsO2a1iLWoZy44QQQgghno/Oid3q1avx8vJi27ZtivJ27drh4+PD4MGDWb16\ntSR2Oih6GtZdrq8TQgghhB7o/LiT27dva948UZSBgQGvvfYat2/f1ltgNVlx74gVQgghhHheOid2\nHTp0IDk5ucT5Fy5coEOHDnoJqqY7fUt5KtbdSkbshBBCCPH8dD4VGxkZiZ+fH02aNGHs2LGYm5sD\nkJ2dzdq1a9m9ezdbt24tt0BrirTsfH6//78bJ0yMoLVKbpwQQgghxPMrMaPw8PDQKjMwMGDOnDnM\nnTsXGxsbDAwMSE9Pp7CwEBsbG8aOHUtCQkK5BlzdxacrR+s6WRlTW26cEEIIIYQelJjYWVlZYVDk\nFVfW1tY4Ozsrypo1a1Y+kdVQx39/oJjubiunYYUQQgihHyUmdrt3767IOP4yio7YdbM1qaRIhBBC\nCFHT6HzzhHh+v98r4GJWvmba0AC62siInRBCCCH046mu2s/Ly2Pjxo3s37+f1NRU4NErxV5++WVG\njRpF7dry2I7S7P8tVzHd0bI2FsaSWwshhBBCP3RO7DIzM/Hx8eHs2bPY2NjQvHlzAE6fPs2BAwfY\nuHEjO3fuRKVSlVuw1d3eVGVi93IT00qKRAghhBA1kc7DRXPmzOHChQtERUVx4cIF9u7dy969e0lO\nTmbVqlVcuHCBjz/+uDxjrdZy8go5dEN548RLktgJIYQQQo90Tuz27NnDuHHjGDFiBIaG/1vMwMCA\nYcOGERwcLDdclGJXSi738tWaaQczI9o3lFPXQgghhNAfnRO7rKysUh9t0qxZM7KysvQSVE305aUc\nxfSQ5nW0HicjhBBCCPE8dE7smjdvzp49e1Cr1Vrz1Go1u3fv1lx3J5SSbj3kyO/Kx5wMdzGrpGiE\nEEIIUVPpnNgFBwcTFxeHn58f+/fv5+rVq1y9epV9+/bh5+fH4cOHCQkJKc9Yq635SXcV091sjWmt\nktOwQgghhNAvne+KDQoKIiMjg0WLFhEXF6cpV6vVGBsb8+GHHzJmzJhyCLF623ntPt+lKe+G/cCt\nXiVFI4QQQoia7KmeYzdlyhSCgoKIi4sjLS0NtVqNo6MjXl5eNGzYsLxi1MnatWtZvnw56enptG7d\nmoiICLp3716pMcX9v1zePvKnouwF69r0d5C3TQghhBBC/54qsQOwtLTEz8+vPGJ5Ztu2bSMsLIxP\nP/0UT09P1q5di7+/PwkJCTRp0qRCYzn/Zx6zk43JSL7Jf2/lKebVMoAl3RvITRNCCCGEKBc6X2N3\n8eJFrceZHDt2DF9fX/r168dnn32m9+B0FRUVxYgRIxg9ejStWrUiMjISW1tb1q9fX+GxZD0sZP+t\nWlpJHUCER33aySNOhBBCCFFODDIzM7Vvcy2Gv78/BgYGfP311wBcv34dDw8PTExMsLa25pdffmHl\nypWMGDGiXAMu6uHDh9jZ2bFu3ToGDRqkKQ8NDeX8+fPs2bNHa5lLly6VWzyXcwwYfqqOoswANW83\nzWNMk/wSlhJCCCGEKJuLi0up83U+FXv69GkmTJigmd6yZQuFhYUcPXoUOzs7hg8fztq1ays8scvI\nyKCgoABra2tFubW1NTdv3ix2mbI65XnUyc6HU+maaU8bY6Z3sqCX3V/nurpLly6Vax9XddJ+ab+0\nX9r/VyXtr/z265zYZWVlYWlpqZk+cOAAvXr1ws7ODoCXXnqJmTNn6j9CHRW9bk2tVlfKtWzWdYz4\nuOUDWjva096yNnZmRhUegxBCCCH+mnS+xs7a2prU1FQAMjMzOXHiBF5eXpr5Dx48KGnRcmVpaYmR\nkZHW6NytW7e0RvEqgomRAd42BQxoYipJnRBCCCEqlM4jdl5eXqxZswYLCwuOHj0KwCuvvKKZn5yc\njIODg/4jLIOxsTHu7u7ExsYqrrGLjY3Fx8enwuMRQgghhKgsOid2M2fO5PLly3z00UcYGxvz8ccf\n4+joCEBubi47duxg6NCh5RZoaSZMmMD48ePp3LkzHh4erF+/nt9//53AwMBKiUcIIYQQojLonNhZ\nW1uzd+9e7ty5g6mpKcbGxpp5arWamJgYGjduXC5BlsXX15fbt28TGRlJeno6bdq04euvv9YknkII\nIYQQfwVP/YBiCwsLrbI6derQvn17vQT0rIKDgwkODq7UGIQQQgghKlOpid1///vfp15h586dnzkY\nIYQQQgjx7EpN7Pr376/zI0MeP17k9u3beglMCCGEEEI8nVITu6ioqIqKQwghhBBCPKdSE7uKfouE\nEEIIIYR4djo/oFgIIYQQQlRtktgJIYQQQtQQktgJIYQQQtQQBpmZmerKDkIIIYQQQjw/GbETQggh\nhKghJLETQgghhKghJLETQgghhKghJLETQgghhKghJLETQgghhKghJLETQgghhKghJLHTo7Vr1+Lm\n5oatrS0vvvgix48fr+yQntvixYvx8vKiSZMmtGjRgoCAAM6fP6+o89Zbb6FSqRT/+vfvr6jz4MED\npkyZQvPmzbG3t2fYsGFcv369IpvyzCIiIrTa17JlS818tVpNREQErVu3plGjRrz66qtcuHBBsY7M\nzExCQkJwdHTE0dGRkJAQMjMzK7opz6R9+/Za7VepVAwdOhQou39Atz6qKo4dO8awYcNo06YNKpWK\nL7/8UjFfX9v73LlzvPLKKzRq1Ig2bdqwYMEC1OrKf/pUae3Py8tj1qxZdO/eHXt7e1q1akVwcDBp\naWmKdbz66qtafxNBQUGKOlX1O1HW9tfX/i4tLY2AgADs7e1p3rw5U6dO5eHDh+XevrKU1f7i9gUq\nlYrQ0FBNnep8TNDlmFfV9wGS2OnJtm3bCAsLY/LkyRw+fJiuXbvi7++vtcOrbo4ePcrYsWPZt28f\nMTEx1KpVi0GDBvHnn38q6vXp04eLFy9q/n3zzTeK+eHh4ezatYt169axZ88e7t69S0BAAAUFBRXZ\nnGfm4uKiaN+TSfuyZcuIiopiwYIFHDx4EGtrawYPHszdu3c1dYKDgzlz5gzffPMN0dHRnDlzhvHj\nx1dGU55abGysou2HDh3CwMCAQYMGaeqU1j+gWx9VFTk5ObRt25b58+dTp04drfn62N537txh8ODB\n2NjYcPDgQebPn8+KFStYuXJlhbSxNKW1/969e5w+fZrQ0FAOHTrEpk2buH79OkOGDCE/P19Rd+TI\nkYq/iSVLlijmV9XvRFnbH55/f1dQUEBAQADZ2dns2bOHdevWERMTw/Tp08u9fWUpq/1PtvvixYt8\n9dVXAIr9AVTfY4Iux7yqvg+QBxTrSb9+/XB1dWX58uWask6dOjFw4EBmzZpViZHpV3Z2No6Ojnz5\n5Zd4e3sDj36d3b59my1bthS7TFZWFs7OzkRFRWlGeX777Tfat29PdHQ0/fr1q7D4n0VERAQxMTHE\nx8drzVOr1bRu3Zpx48ZpfrHev38fFxcX5s6dS2BgIBcvXsTDw4PvvvsOT09PAOLj4/H29uann37C\nxcWlQtvzvBYtWsTy5ctJTk7GzMys1P4B3fqoqnJwcGDhwoWMHDkS0N/2XrduHbNnz+aXX37RHDwj\nIyNZv34958+fx8DAoHIaXETR9hcnOTkZT09Pjh07hqurK/BoxK5t27ZERkYWu0x1+U4U13597O8O\nHDjA0KFDOXv2LI0bNwZgy5YtvPfee1y6dAkLC4vyb5wOdNn+7733HsePH+fEiROaspp0TCh6zKsO\n+wAZsdODhw8fkpSURN++fRXlffv25ccff6ykqMpHdnY2hYWFqFQqRXl8fDzOzs507tyZ9957jz/+\n+EMzLykpiby8PEX/NG7cmFatWlWb/rl27Rpt2rTBzc2NoKAgrl27BkBKSgrp6emKttWpU4fu3btr\n2paYmIi5uTkeHh6aOp6entStW7fatP8xtVrNF198QUBAAGZmZprykvoHdOuj6kJf2zsxMZFu3bop\nRkT69evHjRs3SElJqaDW6MfjUYqi+4StW7fSvHlzPD09mTFjhmI0o7p/J553f5eYmEirVq00SR08\n2v4PHjwgKSmp4hrynLKzs9m2bRujR4/WmldTjglFj3nVYR9Q67mWFgBkZGRQUFCAtbW1otza2pqb\nN29WUlTlIywsjPbt29O1a1dNWf/+/Xn99ddp2rQpqampzJs3Dx8fH+Li4jAxMeHmzZsYGRlhaWmp\nWFd16Z8uXbrw2Wef4eLiwq1bt4iMjGTAgAEkJCSQnp4OUOy2v3HjBgA3b97E0tJS8QvMwMAAKyur\natH+J8XGxpKSksKoUaM0ZaX1T8OGDXXqo+pCX9v75s2b2Nvba63j8TwnJ6fyaoJePXz4kBkzZvDy\nyy/j4OCgKff396dJkyY0atSI5ORk5syZw88//8yOHTuA6v2d0Mf+7ubNm1p/Q5aWlhgZGVX59j8p\nOjqaBw8eMHz4cEV5TTomFD3mVYd9gCR2elR06FStVleZUyr68OGHH5KQkMB3332HkZGRptzPz0/z\nf1dXV9zd3Wnfvj379u3Dx8enxPVVl/7529/+ppju0qUL7u7ubNq0iRdeeAEoe9sX187q0v4nbdy4\nkU6dOuHm5qYpK61/3nnnHU15Tfp+6GN7F7eOkpativLz8wkJCSErK4vNmzcr5o0ZM0bzf1dXV5yc\nnOjXrx9JSUm4u7sD1fc7oa/9XUntrOrtf9LGjRt59dVXsbKyUpTXlGNCScc8qNr7ADkVqwcl/dK6\ndeuWVlZfXYWHh7N161ZiYmLK/CVhZ2eHvb09V69eBcDGxoaCggIyMjIU9apr/5ibm9O6dWuuXr2K\nra0tQKnb3sbGhlu3binudlKr1WRkZFSr9v/xxx/s2bOn2NMuT3qyfwCd+qi60Nf2trGxKXYdoD0S\nUBXl5+czduxYzp07x86dO2nYsGGp9Tt27IiRkZFin1ATvhPwbPu74rZ/SWd+qqozZ85w6tSpMvcH\nUD2PCSUd86rDPkASOz0wNjbG3d2d2NhYRXlsbKziHHt1NW3aNKKjo4mJidF6jEVxMjIyuHHjhuYL\n4O7uTu3atRX9c/36dc0FptVNbm4uly5dwtbWlqZNm2Jra6toW25uLvHx8Zq2de3alezsbBITEzV1\nEhMTycnJqVbt37RpEyYmJvj6+pZa78n+AXTqo+pCX9u7a9euxMfHk5ubq6kTGxuLnZ0dTZs2raDW\nPJu8vDwCAwM5d+4cu3bt0mzn0pw7d46CggJN3ZrynYBn29917dqVixcvKh7vERsbi4mJiWZEs6rb\nuHEjjo6O9OnTp8y61e2YUNoxrzrsA4zCwsJmP9caBAD16tUjIiKCRo0aYWpqSmRkJMePH2flypXU\nr1+/ssN7ZqGhoXz11Vds2LCBxo0bk5OTQ05ODvAooc3Ozubjjz/G3Nyc/Px8zp49y7vvvktBQQGR\nkZGYmJhgamrK77//zr/+9S/atWtHVlYWkyZNwsLCgjlz5mBoWLV/X8yYMQNjY2MKCwu5fPkyU6ZM\n4erVqyxZsgSVSkVBQQFLlizB2dmZgoICpk+fTnp6OkuXLsXExAQrKytOnDhBdHQ0bm5uXL9+nUmT\nJtGpU6cq8XgHXajVaiZMmMBLL72k9ViD0vqnfv36GBgYlNlHVUl2djbJycmkp6fzxRdf0LZtWyws\nLHj48CH169fXy/Zu0aIFn3/+OWfPnsXFxYX4+HhmzpzJ+++/X+kHttLaX7duXUaPHs3Jkyf597//\nTb169TT7BCMjI2rXrs2vv/7KmjVrqFu3Lg8fPiQxMZH3338fBwcHZsyYgaGhYZX+TpTWfiMjI73s\n75ycnNi1axcHDx7E1dWV5ORkQkND8ff35/XXX6+y7X98LLt37x5vv/02ISEh9OjRQ2v56nxMKOuY\np8v+rLL3AfK4Ez1au3Yty5YtIz09nTZt2vCPf/xD64++uil6p9tj06ZNIzw8nPv37zNy5EjOnDlD\nVlYWtra29OrVi+nTpyvu+MrNzeWjjz4iOjqa3NxcevfuzaeffqqoU1UFBQVx/PhxMjIysLKyokuX\nLkyfPp3WrVsDj5Ke+fPns2HDBjIzM+ncuTOLFi2ibdu2mnX8+eefTJs2jb179wLg7e3NwoULS+zf\nqubw4cP4+Pjwww8/0LlzZ8W8svoHdOujquLIkSPFHlyHDx/OqlWr9La9z507R2hoKCdPnkSlUhEY\nGMi0adMq/Rqj0tofFhZGhw4dil0uKiqKkSNH8ttvvxESEsKFCxfIycnBwcGBAQMGEBYWRoMGDTT1\nq+p3orT2L168WG/7u7S0NEJDQzl8+DCmpqYMGTKEefPmVfoPnbL+/gH+85//MHHiRH7++Wfs7OwU\n9ar7MaGsYx7ob59fXvsASeyEEEIIIWqIqn0OTAghhBBC6EwSOyGEEEKIGkISOyGEEEKIGkISOyGE\nEEKIGkISOyGEEEKIGkISOyGEEEKIGkISOyGEEArvvvsugwcPfqZlT506hZWVFVeuXNFzVEIIXUhi\nJ4R4ZiqVSqd/X375ZWWHWiUcO3aMiIgIsrOzKzuUEl25coXNmzczadIkTVlubi4qlYqwsDCt+kuX\nLkWlUhESEkJBQQEdO3akV69eREREVGTYQoj/r1ZlByCEqL5Wr16tmN6wYQMnTpxg5cqVivLKfk1W\nVXH8+HEWLFhAUFAQ5ubmlR1OsaKionBycqJ3795l1l22bBmzZ89m6NChrFq1CiMjIwACAwMZM2YM\ns2fPrvQ3CQjxVyOJnRDimQUEBCim4+LiOHnypFZ5TXXv3j3MzMwqOwxycnKoW7fuc68nNzeX6Oho\nQkJCyqy7YsUKZs2apZXUAQwYMIC6deuyefNmpkyZ8txxCSF0J6dihRAVprCwkJUrV+Lp6YmtrS0u\nLi5MnDiRzMxMRb2WLVsyYsQIYmNj6d27N40aNaJHjx4cP34cgB07dmjW4eXlxdmzZxXLBwUF4ejo\nyK+//oqvry/29va0atWKefPmUVBQoBXXpk2bePHFF2nUqBFOTk4EBQWRlpamqNO/f3969uzJqVOn\n8Pb2xs7OjunTpwOP3qU7evRo2rVrh42NDW3atOGDDz4gKytLs/zs2bP55JNPAGjVqpXmNPVPP/2k\nOdW5ZMkSrdhatmypOC26fv16VCoV8fHxTJ48GWdnZ1q0aKGZn5mZydSpU3F1dcXGxoaOHTuyePFi\nCgsLy9w+R44c4c6dO/Tt27fUeitWrOCjjz4qNqkDMDU1xcPDg2+//bbMzxRC6JeM2AkhKsw777zD\nN998w8iRIxk/fjypqamsWbOGpKQkDhw4gLGxsabu5cuXefPNNwkKCiIgIIDly5czbNgwIiMjmTdv\nHmPHjqWwsJAlS5YQFBREYmKi4uXZ+fn5+Pr60rVrV+bMmUNcXByLFi0iJydHcf1XREQECxcuxM/P\nj7///e9kZGSwZs0avL29OXr0qOKl3RkZGQwZMgRfX18CAgJo2LAhAFu3buXevXsEBgZiaWnJ2bNn\n+eKLL7h48SK7d+8GwNfXl19//ZWdO3cSGRmJhYUFAM2bN3+mvpw0aRKWlpZMnTqVu3fvApCdnc0r\nr7xCeno6gYGBODg4kJiYyNy5c7lx4waRkZGlrvPHH3/E0NAQNze3EutERUWVmtQ95u7uzpIlS8jO\nzq6yp52FqIkksRNCVIhDhw6xadMmPv/8c8Udl3369GHgwIFs3bqV4cOHa8p/+eUXfvjhBzp37gxA\ns2bNGDFiBJMmTeKnn37CwcEBADMzM8LCwkhISKBbt26a5e/du8eAAQNYsGABAOPGjSMwMJDVq1fz\n1ltv4ejoyJUrV4iMjGT27NlMnDhRs+zAgQPp2bMnq1evZtq0aZryGzdusHTpUsaMGaNoW0REhNYp\nWXd3d959911OnTpFx44dcXNzo127duzcuRMfHx9sbW01dXNzc5+6Pxs0aEBMTIwisVq2bBmpqakc\nOXKEZs2aAWgSvMWLFzNhwgScnJxKXOelS5ewsrIqMRHbvXs3aWlpZSZ1AE5OThQUFHD58mXc3d2f\nun1CiGcjp2KFEBVix44dqFQqevfuTUZGhuZfu3btsLCw4MiRI4r6rq6umqQOoEuXLsCjRPBxUvdk\n+bVr17Q+880331RMh4SEUFhYyPfffw9ATEwMarWagQMHKmKysrLCxcVFK6Y6deowcuRIrc95nNSp\n1Wru3LlDRkaG5oaRpKQknfrnaY0ZM0YrsdqxYwc9evTAwsJC0Z4+ffpQWFioOZVdktu3bytGKIv6\n448/gEdJW2lJHaBZz+3bt3VpjhBCT2TETghRIS5fvkxmZqbierAn3bp1SzFd9G7Kx6cun0zqniwv\nep1erVq1aNq0qaLM2dkZQHP93OXLl1Gr1SWOKLVu3Vox7eDgQO3atbXqXbt2jVmzZvHDDz9oPcrk\nzp07xa77eT0ekXtMrVZz9epVLl26pHMfF0etVpc4b9SoUaSmprJw4UJUKhVvv/12met58vS4EKL8\nSWInhKgQhYWFNGrUiH/+85/Fzn98vdpjJY0IlVReWkJSUp3CwkKMjIyIjo4uNgEpenrV1NRUq05e\nXh6DBg0iJyeH0NBQXFxcqFu3Lvfv32f48OE63bRQWvJT3M0excWiVqtRq9X079+fd955p9hlBM9+\n3wAABCdJREFUyrqez9LSkuTk5BLnGxkZsWHDBoYMGcL06dOxsLDgjTfeKLbu40S76HYVQpQvSeyE\nEBWiWbNmJCYm0q1bN0xMTMr98/Lz80lJSVGMbD1+G0KTJk00MRUUFNCsWbNSrz0rTVJSEteuXWP9\n+vX4+vpqys+dO6dVt6QEzsTEBDMzM8VdtPDoMSYZGRk6xWFoaIijoyM5OTn06dNH9wY8oWXLluzc\nubPUGx7q1KnD5s2b8fHxYeLEidSrV4+BAwdq1UtJScHIyKjE0UMhRPmQa+yEEBXC19eXvLw8Fi1a\npDUvPz9f61SqPhQdHVyzZg2Ghob069cPgEGDBmFoaMj8+fO1RvPUarVO14c9HkEsunzRhzTD/0YA\ni2urk5OT1jVw69atK/Pzn+Tr60t8fDyHDh3SmpeVlUVeXl6py3t6eqJWqzl9+nSp9SwsLNi6dSvO\nzs6MGzeOgwcPatU5ffo0rq6uckesEBVMRuyEEBWib9++jBo1isjISM6cOcOLL75I7dq1uXLlCjEx\nMcydOxc/Pz+9fZ6ZmRn79+8nMzOTF154gdjYWHbv3s348eM11961bNmS8PBwPvnkE1JSUvD29sbc\n3JyUlBS+/fZb3njjDcUz5IrTtm1bHB0dmTZtGteuXcPCwoJ9+/aRnp6uVbdjx44AzJw5k8GDB1O7\ndm28vLxo2LAho0aNIjw8nMDAQHr16kVSUhLHjh3TXEOoiw8++ID9+/czZMgQRowYQYcOHcjJyeH8\n+fPExMRw8uRJxd24RXXv3p369esTFxdHjx49Sv0sS0tLtm/fzssvv8wbb7zB9u3bNTeM5ObmkpCQ\nUOIpYSFE+ZHETghRYZYvX06nTp3YsGEDc+fOpVatWjRp0gR/f3/Fo0r0oVatWmzbto3Jkyczc+ZM\nzM3NmTx5MuHh4Yp6U6ZMoWXLlqxatUrzaBQHBwf69u3La6+9VubnmJqasmXLFsLDw1m6dCm1atVi\nwIABLFu2DFdXV0Xdbt26ER4ezsaNGzlw4ACFhYUcOHCAhg0bEhISwm+//cbmzZvZt28fvXr1Yvv2\n7fTv31/nNpubm7N3714WL17Mjh072Lx5M/Xq1cPZ2ZmwsDAaNGhQZlv8/f3Zvn275uHLpbG3t2fH\njh14e3szdOhQdu3ahZubGwcOHCA7O1vx+BohRMUwyMzMLPuKYyGEqEaCgoL4/vvvSU1NrexQqp0r\nV67g4eHBtm3bdHpfbHEGDx5MgwYNWL9+vZ6jE0KURUbshBBCaLRo0YIRI0awePHiZ0rskpKSOHLk\nCAkJCeUQnRCiLJLYCSGEUFi+fPkzL+vu7q7T8/KEEOVD7ooVQgghhKgh5Bo7IYQQQogaQkbshBBC\nCCFqCEnshBBCCCFqCEnshBBCCCFqCEnshBBCCCFqCEnshBBCCCFqiP8DTHLtSb6WIWcAAAAASUVO\nRK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1126ba2b0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Note: 4 moles of atoms per formula unit (Fe3C1). That's why we multiply times 4\n",
    "fig = plt.figure(figsize=(9,6))\n",
    "fig.gca().set_xlabel('Temperature (K)')\n",
    "fig.gca().set_ylabel('Isobaric Heat Capacity (J/mol-formula-K)')\n",
    "fig.gca().plot(result['T'], np.squeeze(4.0 * result['heat_capacity']))\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}
